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We  present  numerical  calculations  of  a  four-point  dynamic  susceptibility,  xddi  f°r  the  Kob- 
Andersen  Lennard-Jones  mixture  as  a  function  of  temperature  T  and  density  p.  Over  a  relevant 
range  of  T  and  p,  the  full  t-dependence  of  X4(t)  and  thus  the  maximum  in  x-i(t),  which  is  proportional 
to  the  dynamic  correlation  volume,  are  invariant  for  state  points  for  which  the  scaling  variable  p7/T 
is  constant.  The  value  of  the  material  constant  7  is  the  same  as  that  which  superposes  the  relaxation 
time,  t,  of  the  system  versus  p7/T.  Thus,  the  dynamic  correlation  volume  is  a  unique  function  of 
r  for  any  thermodynamic  condition  in  the  regime  where  density  scaling  holds.  Finally,  we  examine 
the  conditions  under  which  the  density  scaling  properties  are  related  to  the  existence  of  strong 
correlations  between  pressure  and  energy  fluctuations. 
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During  the  last  decade  extensive  evidence  has  accu¬ 
mulated  that  the  dynamics  of  molecules  in  supercooled, 
highly  viscous  liquids  is  inherently  heterogeneous  jj],[§,  3]; 
that  is,  spatial  variations  in  mobility  persist  for  times 
commensurate  with  the  structural  relaxation  time  r. 
Thus,  there  is  growing  interest  in  characterizing  the  dy¬ 
namic  heterogeneities  to  better  understand  the  phenom¬ 
ena  associated  with  the  many-body  dynamics  of  vit- 
ryifing  liquids.  These  phenomena  include  rotational- 
translational  decoupling  HI,  the  dynamic  cross-over 
[d!,  [7J,  enhanced  mobility  under  confinement  non- 
exponentiality  of  the  relaxation  [9|,  and  most  promi¬ 
nently,  the  slowing  down  of  the  translational  and  rota¬ 
tional  dynamics  upon  approach  to  the  glass  transition 
temperature  Tg.  To  explain  these  features,  theories  of  the 
glass  transition  as  diverse  as  the  classical  Adam-Gibbs 
entropy  0]  and  Cohen- Grest  free  volumes  M  models, 
as  well  as  more  modern  approaches  0,0,  invoke  dy¬ 
namic  heterogeneity  having  a  length  scale  that  grows  in 
concert  with  r  upon  drawing  near  to  Tg. 

A  proper  description  of  the  dynamic  heterogeneities  re¬ 
quires  multi-point  dynamic  susceptibilites,  which  reflect 
correlations  in  the  spatial  variation  of  the  dynamics.  A 
four-point  dynamic  susceptibility  x*{t)  can  be  calculated 
as  the  variance  of  the  self-intermediate  scattering  func¬ 
tion  Fs(k,  t) 

X4  (t)=N[(f*(k,t))-F*(k,t)]  (1) 

where  fs(k,t)  is  the  instantaneous  value,  such  that 
( fs(k ,  t))  =  Fs(k,  t).  Xi(t)  quantifies  the  amplitude  of  the 
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fluctuations  associated  with  /s(/c,  t)  and  has  a  maximum, 
xTax  =  X4  (tmwrb  proportional  to  the  dynamic  correla¬ 
tion  volume  [14  HE].  exhibits  various  regimes  m, 

but  most  interesting  for  study  of  the  glass  transition  is 
the  behavior  around  tmax  ~  r  0].  Recent  numerical  and 
simulation  works  have  shown  that  the  dynamic  correla¬ 
tion  volume  grows  upon  cooling  mim.  However, 
very  little  is  known  about  the  combined  temperature  and 
density  dependences  of  %4(t).  In  particular,  unexplored 
is  the  possibility  of  a  description  of  the  p  and  T  depen¬ 
dences  of  Xi{t)  iu  terms  of  the  scaling  property  estab¬ 
lished  for  the  structural  relaxation  time  :20.  l2ll  i22| 

t  =  W/T)  (2) 

where  T\  is  a  function  and  7  a  material  constant.  Experi¬ 
ments  have  shown  that  Eq.  ©,  with  similar  relations  for 
the  diffusion  constant  and  viscosity,  applies  universally 
to  organic,  non-associated  liquids  [23| ,  with  a  range  of 
validity  extending  from  the  high  temperature  Arrhenius 
regime  down  to  Tg  0] .  We  anticipate  a  similar  relation 
for  tmax  and  examine  the  possibility  that  this  scaling 
property  extends  to  the  dynamic  correlation  volume 

xTx=wm  (3) 

To  cast  this  work  in  more  general  terms,  we  assess 
whether  for  a  prototypical  model  glass-former,  the  den¬ 
sity  scaling  properties  hold  for  the  full  time  dependence 
of  both  Fs(k,t)  and  X4(i),  using  the  same  scaling  expo¬ 
nent  7.  The  connection  between  the  dynamic  properties 
of  viscous  liquids  and  7  is  intriguing  because  molecular 
dynamics  (MD)  simulations  have  shown  this  parameter 
to  be  a  measure  of  the  correlation  between  fluctuations 
in  the  potential  energy  and  the  virial  0,  [2(|  0]  ,  pos¬ 
sibly  reflecting  a  hidden  scale  invariance  in  viscous  liq¬ 
uids  0.  7  is  also  related  to  the  steepness  of  the  effec¬ 
tive  repulsive  potential  in  the  range  of  closest-approach 
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FIG.  1:  Reduced  relaxation  times  t*  =  r{p1^3T1^2)  as  a  func¬ 
tion  of  p1  IT  with  7  =  5.1  for  all  studied  densities.  Inset: 
self  intermediate  scattering  functions  as  a  function  of  reduced 
time  t*  =  t(p1^3T1^2)  for  state  points  at  which  p7/T  =  5.07: 
T  =  0.402  (at  p  =  1.15),  T  =  0.50  (at  p  =  1.20),  T  =  616  (at 
p  =  1.30),  and  T  =  0.912  (at  p  =  1.350).  A  constant  reduced 
wave-vector  k*  =  k(p1^3)  =  7.44  is  considered. 


between  particles  [29|,  [30|,  |31j ,  although  this  connection 
may  seem  tenuous  given  the  limitations  of  a  two-body  po¬ 
tential  in  describing  interactions  in  real  liquids  and  the 
possible  non-trivial  role  of  attractive  forces  (31,  32i.  '33  . 
We  further  analyze  this  herein  by  evaluating  pressure- 
energy  correlations  for  the  model  studied  in  1331] .  which 
did  not  conform  to  density  scaling  over  the  relevant  den¬ 
sity  regime. 

In  this  work  we  consider  the  well-studied  Kob- 
Andersen  (KA)  binary  mixture  [1H  as  a  model  glass¬ 
forming  liquid.  It  consists  of  1000  particles  in  a  cubic 
box  with  periodic  boundary  conditions.  Particles  inter¬ 
act  through  the  Lennard-Jones  (LJ)  potential 


uap{r)  =  4ea/3 


(4) 


where  a,  (3  =  1,2  are  species  indices.  The  values  of  the 
parameters  in  Eq.  ((4j)  can  be  found  in  the  original  paper 
[341].  In  the  following  we  use  reduced  LJ  units,  assuming 
(Jii,  en,  and  ^/micr^/en  (where  mi  is  mass)  as  units  of 
distance,  energy,  and  time,  respectively.  We  performed 
MD  simulations  in  the  NVT  ensemble  using  the  Nose- 
Poincare  thermostat  [35[  with  a  mass  parameter  Q= 5.0. 
We  considered  five  isochoric  paths  in  the  density  range 
1.150  <  p  <  1.350.  For  each  state  point  we  averaged  the 
dynamic  properties  over  20  independent  realizations  of 
the  system. 

Following  previous  work  and  to  emphasize  the  con¬ 
nection  with  the  exact  density  scaling  relations  observed 
for  inverse  power  law  (IPL)  potentials,  we  analyze  “re¬ 
duced”  quantities  (indicated  by  stars)  using  p-1/3  and 
T1/2  as  reduction  parameters  for  distances  and  veloc¬ 
ities,  respectively.  We  begin  our  investigation  of  the 


FIG.  2:  Maximum  of  four  point  dynamic  susceptibility  as  a 
function  of  p7 /T  with  7  =  5.1  for  all  studied  densities.  Inset: 
four  point  dynamic  susceptibility  as  a  function  of  reduced  time 
t*  for  state  points  at  which  p7/T  =  5.07  (same  state  points 
as  inset  of  Fig.  [lj . 


density  scaling  properties  of  the  KA  model  by  calcu¬ 
lating  Fs{k,t)  =  (4/^0  SyLi(exP (*k  •  [rj(t)  —  (0)] } ) 
for  the  various  densities  at  a  fixed  reduced  wave-vector 
k*  =  ^(p1/3)  =  7.44,  which  matches  the  position  of  the 
first  peak  in  the  static  structure  factor  ( k  =  7.0)  for  the 
well-studied  density  p  =  1.2.  The  reduced  relaxation 
times  r*,  defined  as  the  time  for  Fs{k*,t*)  to  decay  by  a 
factor  of  e,  are  shown  in  Fig.  [T]  versus  the  scaling  variable 
p7/T.  The  material  constant  7  =  5.1  ±  0.1  provides  the 
optimal  collapse  of  r*  onto  a  single  curve;  this  value  is  in 
accord  with  the  scaling  behavior  found  previously  for  the 
diffusion  coefficient  of  this  system  t30| .  We  note  that  the 
pressures  attained  at  the  lowest  temperatures  range  from 
1  to  20  reduced  LJ  units,  depending  on  p.  Assuming  Ar¬ 
gon  units,  this  corresponds  to  0.04  to  0.8  GPa,  which  is  a 
significant  and  experimentally  accessible  pressure  range. 

It  has  been  demonstrated  experimentally  that  the 
shape  of  the  frequency-dependent  linear  response  func¬ 
tion  is  a  function  of  r  and  thus  of  p7/T  [36|,  )33|,  a  result 
consistent  with  the  existence  of  “isomorphic”  points  in 
liquid  state  diagrams  (38|.  We  calculated  the  interme¬ 
diate  scattering  functions  for  each  of  the  five  densities 
at  the  respective  temperatures  corresponding  to  a  fixed 
value  of  p7/T  =  5.07.  In  the  inset  of  Fig.  |T]  these  cor¬ 
relation  functions  are  plotted  as  a  function  of  reduced 
time,  and,  as  found  previously  over  a  more  limited  den¬ 
sity  range  (38),  Fs(k*,t*)  has  essentially  the  same  shape 
for  state  points  for  which  p7/T  is  constant.  Thus,  not 
only  do  the  relaxation  times  superpose  as  a  function  of 
p7 /T,  but  the  entire  t-dependence  of  the  correlation  func¬ 
tions  is  invariant  for  isomorphic  state  points. 

We  now  examine  the  correlations  in  the  dynamic  fluc¬ 
tuations  to  assess  explicitly  whether  density  scaling  can 
be  extended  to  high-order  correlation  functions,  as  envis¬ 
aged  in  [js|.  To  do  this  we  calculate  the  four  point  dy- 
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FIG.  3:  Maximum  of  four  point  dynamic  susceptibility  as  a 
function  of  reduced  relaxation  times  r*  for  all  studied  densi¬ 
ties  p.  Inset:  reduced  time  at  which  \4 (i*)  is  maximum  as  a 
function  of  t*. 


namic  susceptibility  [Eq.  (JT|)]  associated  with  the  complex 
instantaneous  value  of  fs(k,t)  at  the  same  fixed  reduced 
wave-vector  k*  considered  above.  The  use  of  the  com¬ 
plex  self  intermediate  scattering  function  (rather  than 
the  real  part  used  in  previous  studies)  removes  the  finite, 
long-time  limit  from  x<t{t)  without  altering  the  general 
features  of  the  correlation  function.  In  Fig.  [2]  the  maxi¬ 
mum  of  X4(t)  is  plotted  as  a  function  of  p1  /T  using  the 
aforementioned  7  =  5.1.  To  provide  an  estimate  of  the 
statistical  uncertainties,  we  include  errors  bars  that  rep¬ 
resent  two  standard  deviations  on  the  average  over  sys¬ 
tem  realizations  for  selected  states.  As  seen  in  Fig.  [2] 
density  scaling  applies  to  xTax  within  the  estimated  er¬ 
ror,  using  the  same  value  for  the  scaling  exponent  that 
superposes  the  relaxation  times.  In  the  inset  of  Fig.  [2] 
we  show  the  full  t-dependence  of  X-rW  f°r  state  points 
at  which  p1  /T  =  5.07.  As  is  the  case  for  the  “average” 
intermediate  scattering  functions,  within  the  estimated 
error  bars  the  X4(t)  fall  011  a  single  curve.  Deviations 
from  the  scaling  are  observed  only  for  the  lowest  investi¬ 
gated  density  (p  =  1.15). 

To  make  contact  with  previous  numerical  and  exper¬ 
imental  work  on  four  point  dynamic  susceptibilities,  we 
show  in  Fig.  [3]  xTax  as  a  function  of  the  reduced  re¬ 
laxation  times,  r* .  Since  xTax 's  proportional  to  the  dy¬ 
namic  correlation  volume,  Fig.  [3] confirms  the  presence  of 
a  steady  but  rather  mild  growth  of  dynamic  correlations 
as  the  structural  relaxation  times  increase  to  the  point  of 
vitrification.  It  also  shows  that,  at  fixed  r*,  the  dynamic 
correlations  are  invariant  to  either  T  or  p,  as  expected 
from  the  density  scaling  of  both  the  average  dynamics 
and  the  dynamic  correlations  for  the  same  value  of  7. 
Recent  experiments  0  have  shown  that  a  similar  result 
holds  for  the  temperature  derivatives  of  the  two  point  dy¬ 
namic  correlation  function,  xt  (t) ,  which  provides  a  lower 
bound  to  X*i(t)  00.  This  correspondence  supports 


FIG.  4:  Slope  of  U-W  correlations  as  function  of  T  at  p  =  1.2 
(circles)  and  p  =  1.3  (squares)  for  the  KA  LJ  model  (open 
symbols)  and  its  WCA  variant  (filled  symbols).  The  dashed 
lines  indicate  the  range  of  the  dynamic  scaling  exponent  7  = 
5.1  ±0.1.  Inset:  Pearson  correlation  coefficient  as  a  function 
of  T  for  p  =  1.2. 


the  validity  of  the  experimentally  accessible  Xr(t)  as  an 
approximation  to  %4 (t)-  In  the  inset  of  Fig.  [3l  we  show 
the  r*-dependence  of  the  reduced  time  associated  with 
the  maximum  in  Xi{t)-  As  expected,  the  two  quantities 
are  essentially  equal. 

Recent  numerical  work  [2Sl.  [2(1  [27j  has  shown  that  the 
dynamic  scaling  exponent  7  can  be  independently  esti¬ 
mated  on  the  basis  on  the  correlation  between  fluctu¬ 
ations  of  two  thermodynamic  quantities,  the  potential 
energy,  U,  and  the  virial,  W  (the  configurational  part 
of  the  pressure).  These  fluctuations,  AU  =  U  —  ( U ) 
and  AW  =  W  —  (W)  are  proportional  for  particles  in¬ 
teracting  with  IPL  potentials  [25|  and  have  been  shown 
to  be  strongly  correlated  (Pearson  correlation  coefficients 
R  >  0.9)  for  various  other  liquids  Jig!].  In  the  latter  cases 
the  slopes,  obtained  from  linear  regression  of  AW  vs. 
AU,  are  equal  within  the  statistical  fluctuations  to  the 
dynamic  scaling  exponent  [22}  ■  Such  results  support  the 
conjecture  that  liquids  display  strong  U-W  correlations 
if  and  only  if  they  comply  with  density  scaling  [25],  the 
inference  being  that  these  properties  have  a  common  ori¬ 
gin  in  the  same  generalized  IPL  approximation  of  the 
interaction  potential  [26j. 

To  further  elaborate  on  these  aspects,  we  show  in  Fig. [4] 
the  AW  vs.  AU  slopes,  T,  obtained  for  the  KA  model, 
as  a  function  of  T  for  two  densities  ( p  =  1.2  and  1.3). 
r(T,  p)  is  always  close  to  the  scaling  exponent  (within  the 
estimated  uncertainty),  as  found  previously  from  simu¬ 
lations  along  isobaric  paths  [27)].  Also  included  in  Fig.  [4] 
are  results  from  additional  simulations  carried  out  using 
the  purely  repulsive,  Weeks-Chandler-Andersen  (WCA) 
variant  of  the  KA  model  [391  l4p| .  In  the  WCA  model, 
the  interaction  parameters  are  unchanged  but  each  of  the 
pair  potentials,  uap(r),  is  shifted  so  that  the  minimum  is 
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zero  and  the  potential  is  truncated  at  this  minimum  [39j . 
Very  recently,  Berthier  and  Tarjus  [3SS]  concluded  from 
simulations  of  this  model  that  density  scaling  surpris¬ 
ingly  requires  the  contribution  of  the  attractive  interac¬ 
tions;  in  the  WCA  model,  in  fact,  the  scaling  of  r  was 
absent  except  at  very  high  densities.  As  seen  from  Fig.  [4] 
and  its  inset,  in  the  WCA  model  there  are  strong  U-W 
correlations  ( R  >  0.9);  however,  the  value  of  F  changes 
significantly  as  the  state  parameters  are  varied.  That  is, 
the  local  scaling  exponents  change  too  much  (for  reasons 
yet  unknown),  causing  a  breakdown  of  the  density  scal¬ 
ing.  Nevertheless,  for  every  state  point  U  and  W  are 
strongly  correlated. 

In  the  light  of  these  results,  the  conjecture  of  Pedersen 
et  al.  [25]  connecting  strong  U-W  correlations  and  den¬ 
sity  scaling  must  be  partly  reformulated.  The  existence 
of  strong  U-W  correlations  is  not  a  sufficient  condition 
for  Eq.  (J2J)  to  apply.  There  is  an  additional  requirement, 
that  the  slopes  of  the  U-W  correlation  must  be  (almost) 
insensitive  to  variations  of  the  state  parameters.  This  is 
consistent  with  previous  simulation  and  experimental  re¬ 
sults  on  local  dynamic  scaling  exponents  [41],  |42j .  Under 
these  conditions,  density  scaling  applies  to  a  very  good 
approximation  to  all  time-dependent  properties,  includ¬ 
ing  high-order  time-dependent  correlation  functions,  as 
conjectured  in  [38].  On  the  other  hand,  further  work  is 
needed  to  understand  under  which  conditions  the  slopes 


of  U-W  correlation  are  insensitive  to  variations  of  p  and 
T,  and  how  this  is  related  to  the  attractive  part  of  the 
interaction  potential. 

From  simulations  over  a  range  of  T  and  p,  the  dynamic 
correlations  in  a  supercooled  LJ  mixture  are  shown  to 
have  a  spatial  extent  that  depends  only  on  the  quantity 
p7/T.  In  consideration  of  the  general  behavior  of  liquids 
obeying  density  scaling,  this  means  the  dynamic  correla¬ 
tion  volume  is  related  to  other  mutually  correlated  prop¬ 
erties:  the  relaxation  time,  the  shape  of  the  relaxation 
function,  U-W  correlations,  and  at  least  approximately 
the  isobaric  fragility.  From  experiments  it  is  also  found 
that  the  relaxation  time  is  constant  both  at  the  onset 
of  non- Arrhenius  behavior  at  high  T  [43]  and  at  the  dy¬ 
namic  crossover  at  ca.  1.2Tg  (dJE].  Since  xTax  depends 
only  on  r,  these  changes  in  the  dynamics  also  occur  at 
a  fixed  (pressure-independent)  correlation  volume.  Al¬ 
though  examining  repulsive  exponents  other  than  the  LJ 
value  of  12  remains  for  future  work,  if  the  equivalence  of 
the  scaling  exponents  for  and  r  is  maintained,  we 

expect  that  the  dynamic  correlation  volume  should  vary 
among  different  materials. 
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